New positive small vacuum region gravitational energy expressions 
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We construct an infinite number of new holonomic quasi-local gravitational energy-momentum den- 
sity pseudotensors with good limits asymptotically and in small regions, both materially and in 
vacuum. For small vacuum regions they are all a positive multiple of the Bel-Robinson tensor and 
consequently have positive energy. 
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The localization of energy for gravitating systems has 
been an important fundamental problem since before 
Einstein finalized his field equations (see [![). In re- 
cent decades there has been some significant progress on 
this outstanding problem. New perspectives and ideas on 
gravitational energy continue to be considered as note- 
worthy (e.g., [HUE], HI)- Recall that it is a well-known 
inevitable consequence of the equivalence principle that 
there is no co variant (reference frame independent) de- 
scription for the gravitational energy-momentum density 
(for a discussion see Ch. 20 in 0]). The coordinate de- 
pendent pseudotensor density description used in earlier 
times has in recent decades largely been replaced by the 
quasi-local perspective (i.e., energy- momentum is to be 
associated with a closed 2-surface; for a comprehensive 
review see Q). However it has been noted that the 
Hamiltonian approach to quasi-local energy-momentum 
includes all the possible pseudotensors simply by tak- 
ing their associated superpotentials as the Hamiltonian 
boundary term. This not only shows the pseudotensors 
to be a special type of quasi-local expression but, more- 
over, reveals their ambiguities to be just the same as 
those of the quasi-local Hamiltonian boundary term. The 
Hamiltonian approach, however, gives these ambiguities 
a clear physical and geometric meaning. Specifically, the 
freedom in the choice of expression is associated (via the 
boundary term in the Hamiltonian variation) with the 
freedom to choose the type of boundary conditions, and 
the choice of quasi-local reference (which determines the 
zero energy or ground state) is effectively just the choice 
of coordinates on the boundary for the holonomic pseu- 
dotensor @, i, El El El . 

Various criteria have been proposed for quasi-local 
quantities (see, e.g., @,EI[), including having the correct 
asymptotic limits at infinity and positivity. Positivity 
has been regarded as a strong condition. It is difficult to 
prove the positivity of the quasi-local energy determined 
by some expression for a general region. Just consid- 



ering a small vacuum region already gives a significant 
requirement, namely that the energy-momentum be pro- 
portional to the Bel-Robinson tensor 0, EH • This will 
guarantee positive energy within a small vacuum region. 

We have had some hope that positivity within a small 
vacuum region would be a quite strong restriction, per- 
haps even selecting a unique expression. We found that 
it is indeed a serious constraint; in particular the re- 
quirement of a positive energy density eliminates all of 
the classical pseudotensors, nevertheless we found that 
it still allows a surpri sing amount of freedom. Specifi- 
cally, we recently [l5|, Efl E3, EH showed that none of 
the classical pseudotensors (Einstein Papapetrou 
2C, [U [12], Landau-Lifshi tz [23l| . Bergman-Thompson 
2|, M0ller [H, Goldberg 0, Weinberg ^) yields a 
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positive energy density within a small vacuum region, 
although certain linear combinations of them do. 

Here we show how to construct an infinite number 
(more precisely an 11 parameter class) of new quasi-local 
pseudotensorial expressions which, in addition to having 
the proper asymptotic and small material limit, have a 
small region energy-momentum density proportional to 
the Bel- Robinson tensor. This guarantees a locally pos- 
itive energy density for small vacuum regions. We will 
present our new results after reviewing some technical 
background. 

Riemann normal coordinates (RNC). As Riemann first 
argued (see, e.g., [28j]), at any preselected point one can 
choose coordinates such that at the point, (i) = 0, 

(ii) the metric coefficients have the standard flat values, 

(iii) the first derivatives of the metric vanish, (iv) and 
the second derivatives have the minimum number (20 for 
n = 4) of independent values. Specifically g a p\o = g a /3, 
9 M 5a/3|o = 0, 39 M „3 a/3 | = -{R aiXr p v + RavPn)\i), where 
R a Piiv is the Riemannian curvature tensor (all indicics 
are holonomic, our conventions, unless otherwise speci- 
fied, follow MTW @), and, in our case, g a p = r\ a fi — 
diag(— 1, +1, +1, +1) is the Minkowsi spacetime met- 
ric. The corresponding Levi-Civita connection values are 

r a /3 7 |o = 0, 39 M r a ;3„|o = —(R a /3v^ + R a v/3n)\o- 

Quadratic curvature basis. It turns out that when 
expanded in RNC the lowest non-vanishing vacuum 
energy-momentum expressions are of the second or- 
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der and are quadratic in the curvature tensor: i„„ 



(R...R.\ 



An investigation [14( of all such pos- 



sible terms (taking into account the Weyl = vacuum Ric- 
mann tensor symmetries) shows that they can be written 

in terms Of Q^aufi '■= R al ibaR a v b p = Qvfijia = Q af j,/3v, 

where all the symmetries have been indicated. The cited 
work defined the three basis combinations X^ap := 

along with the trace tensor T^ap := — \g^Q a a pa- One 
could write all our expansions as linear combinations of 
X, Y, Z, T. However it is more suitable for physical 
purposes to use the Bel-Robinson tensor. 

The Bel-Robinson and some related tensors. The Bel- 
Robinson tensor has many well known remarkable prop- 
erties, see e.g., [HI, H§|- For our considerations we are 
interested in it only in the vacuum, where the Riemann 
tensor reduces to the Weyl tensor. In this case the Bel- 
Robinson tensor is completely symmetric and traceless. 
In vacuum the Bel-Robinson tensor B and two other con- 
venient tensors S and K are given in terms of the afore- 
mentioned basis for quadratic terms by 



B 



Z+3T, S := -2X+2Z-6T, K := Y+9T. (1) 



(Here and below we often suppress obvious indicies.) Our 
leading order non- vanishing vacuum expressions could be 
given as linear combinations of B, S, and K. As can 
be directly verified, these three combinations satisfy the 
second order divergence free condition 



dpi^xh^) = 2x j t 0j a ^ ee 0, 



(2) 



and all such tensors are some linear combination of these 
three. While the tensor 5* has been known for a long time, 
the tensor K (which was recently introduced [la. [l7l. [l8j ) 
also enjoys this divergence free property. However, for 
our purposes here we have found it even more convenient 
to use instead of K the tensor 



V 



-2X + Y + 2Z + 3T ee S + K, 



(3) 



which, in addition to enjoying the vanishing divergence 
property, as we shall see also shares an important positive 
energy property with the Bel- Robinson tensor [16[ . 

All possible superpotentials. One can always express a 
conserved energy-momentum expression in terms of some 

suitable superpotential in the form [H, n nana mm 



(4) 



the anti-symmetry guarantees the conserved current con- 
dition: dpX^v — 0- Properly the superpotential should be 
a density of weight one and the energy-momentum den- 
sity should be a mixed tensor density of the indicated 
type. The physical constraints are such that in the weak 
field, asymptotic, and material limits we want U to be 
linear in T or, equivalently, dg. We consider only su- 
perpotentials with this property. Using the metric and 
connection one can form only three linearly independent 
terms: 
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hl^ aa g a \ \g\h\?g^T a a7 . (5) 



One of these simple choices, 2\g\i g^T^ p u , is the well- 
known M0ller (25[ superpotential, which has an incorrect 
small region matter interior limit [la ]. 

Consider linear combinations of the three terms. In the 
small region matter interior limit the derivative of each 
term gives rise in zeroth order to a certain linear combina- 
tion of IgpiJ^j, and \g\^RS^^. The desired limit to this 
order is 2\g\?G'- L v . (Using Einstein's equation this will 
become 2K\g\^T^ u , the limit required by the equivalence 
principle.) This requirement imposes two constraints on 
the three parameters, hence there is a one parameter set 
which is satisfactory in this limit, namely 



2\g\^ \t^ x K + (3 - 2fc)r 7[ v^ + *4 M r A]7 n 



(0) 



However consideration of obtaining the proper energy- 
momentum asymptotic limit (described in §20.2 in [6|) 
removes this freedom, giving (uniquely up to terms which 
vanish in this limit) a familiar expression, Freud's super- 
potential [13] for the Einstein pseudotensor: 



(Li A 
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= 2 



g \\ fato*\ v + r^v* 1 + <^r A ]-0 , (7) 



which has good limits both at infinity and for small re- 
gions inside of matter. 

Using only the metric and its first derivatives, of the 
classic pseudotensors, one can construct only the Freud- 
Einstein and M0ller superpotential expressions. However 
it should be noted that these expressions really only make 
sense in nearly Cartesian coordinates, hence there actu- 
ally is an implicit background Minkowski structure, a 
reference geometry g,T. If we explicitly introduce this 
structure we can go much further. In this work we will 
use coordinates which are asymptotically flat Minkowski 
for large and small distances. From now on indicies will 
be transvected only with the flat metric g. 

A general superpotential can be expressed as an ex- 
pansion in terms of powers of h := g — g in the form 
U = Uq + Ui(h) + U~2(h 2 ) + • ■ • where Ut is linear in 
r ~ dg. It is easy to see that for the limits of interest we 
need not consider the terms for k > 2. It is well known 
that the Freud superpotential ([7]) gives good values to 
linear order. It has the correct limit both at infinity and 
within matter to zeroth order. Without loss of generality 
we can take U$ to be Freud (alternate choices would just 
shift our expression for U\). Now we want to consider 
expressions of the form 



rrl^A] _ F [A/*](y)(bc) l. ro 



(8) 



where F is some constant tensor made from g a p and 
its inverse with the indicated symmetries. We give 
all possible terms to this order. There are 7 terms 
which include traces: h^T a ^ a , h u ^T^ a a , h a a T^ v , 
h a a 5^Y a X]a , ^ A]c r a ac , S [ ^ a T ab b , h\ L 5 [ ^\, con- 
sequently at second order they give rise to terms pro- 
portional to the Ricci tensor, and hence they vanish in 
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vacuum. They need not be given any more detailed anal- 
ysis here. However there are 6 more terms, 



U A :=2h a ^T a x K, 
U E := 2h bc S^T x \ c , 



U B := 2h b ^T x \ u , 

U D :=2h ba T ab H x \ 

U f := 2h bu T^ b , (9) 



which give rise to the indicated non-vanishing vacuum 
contributions to the energy-momentum density at the 
second order: 



du A ~ 






-1- 2V1 ■ M 


(10) 


dU B ~ 








(11) 


dU c ^ 




^- 


- V\ij^v^ 


(12) 


dU D ~ 






- V]ij^ v , 


(13) 


du E ~ 




3S'4 


2V] ij fl v , 


(14) 


dU F ~ 
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(15) 



(It turns out that two of the expressions produce propor- 
tional results to this order.) 

Any linear combinations of the expressions © , aU A + 
bUB + cUc + dUo + eUE + fUf, may be added to the Freud 
expression. The Freud expression (normalized to give the 
correct asymptotic and small region non- vacuum values) 
leads to a small region vacuum second order contribution 

of sBm 

d x U^K ^ i£[4B - S}^. (16) 

The resultant total energy-momentum expression in 
vacuum to second order then has the form 



where 



(3 = 4-a-b + 2c + 2(d + 2e) + 3f 
1 

s = -1 

v = 2a - b - c+ (d + 2e) 



1 -a-b + 1 -c-l(d + 2e) + 3 -f, 



(17) 

(18) 
(19) 
(20) 



The usual requirement is that one should have a pos- 
itive multiple of the Bel- Robinson tensor [ijj], which 
means that the parameter combinations s and v should 
vanish, and /?, the coefficient of B, should be positive. 
This means, in addition to one inequality, 2 restrictions 
on the 6 parameters controlling the contributions from 
©. Thus there are an infinite number of expressions, 
controlled by 11 parameters (4 free parameters associ- 
ated with these expressions plus 7 additional parameters 
associated with the expressions that explicitly included 
traces), which have positive energy to this order. 

Explicitly we can eliminate and take say a, b, d, e as 
parameters. We then find the two parameter constraints 



3/ = 2 + 3(6 - a) + 2(d + 2e), c = 2a - b + (d + 2e) 
and the inequality needed for positivity is d + 2e > — 1 . 
One simple choice, a = b = c = d= e = and / = 
2/3, yields the one case which had been identified some 
time ago |14| . normalized here to give both the proper 
small region matter interior value as well as the correct 
asymptotic limit. This choice results in a small region 
vacuum energy-momentum density of 



2/e*"" 



> r ? 



3 10 2 



+ 0(x 3 ). 



(21) 



We have a slightly more general result which can conve- 
niently be expressed in another way. Our parameterized 
good superpotential can be put in the form 



U = [t/ FE - (U c + U D )] + (3 



■(Uc + U D ) + U f 



e[U E - 2U D ] 
[U A + 2U C - U F ] + b[U B -U c + U F ] . (22) 



2(U A -U B )--(Uc-U D ) 



Here the first bracket gives the correct material and 
asymptotic limit while its 2nd order vacuum contribution 
vanishes; the second bracket gives a vanishing asymptotic 
and material limit and the desired Bxx/2 in vacuum; the 
third bracket gives Vxx/2 in vacuum, we retain this part 
but have required the Sxx/2 contribution to vanish. The 
remaining 3 terms vanish to these orders, as do the seven 
terms which one could also include built from the afore- 
mentioned explicit trace expressions. 

The associated energy-momentum density for a small 
region, accurate to zeroth order in matter and to second 
order in vacuum, is 



,1 



1, 



The energy-momentum within a small coordinate 
sphere in vacuum is 



P fl = (-E,P) = 



T°„dh 



(24) 



Using 



x a x b d 3 x = -S ab [ r 2 d 3 x = — 5 ab r\ (25) 
3 ./ 3-5 ' V ; 



with a, b, c = 1, 2, 3, we find 
4tt 



2-3-5 



{f3B + vV) a a \ 



(26) 



47T 



= (P + v) — (E ab E ab + H ab H ab , 2e cba E bd H a d ), 

where we have used the remarkable recently discovered 
identity fl6{ . V a a0 v = B a a0 u , along with the traceless 
property, = = B 00 u + B a a0 v , and the well- 

known values of certain components of the Bel- Robinson 
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tensor expressed in terms of the electric and magnetic 
components of the Weyl tensor. From this result it is 
apparent that we are guaranteed to have positive energy 
(in fact a non-spacelike energy-momentum 4-vector) in 
the small sphere limit to second order as long as (3+v > 0. 

Going to higher orders one can expect to find a simi- 
lar situation. There will be a few positivity constraints 
on the terms in the power series expansion, but there 
will also be many additional parameters controlling the 
amplitudes of the numerous possible higher order terms 
in the superpotential. Indeed we can readily see how to 
construct 13 second order in terms just by multi- 
plying the first order ones by h a a . Moreover we can in 
these 13 for any occurrence of /i M „ make the replacement 
hfj, u — > h^-yh? That will generate many more terms. 
Looking to higher orders, we can make similar adjust- 
ments in Ukh k to get many Uk+ih k+1 terms. It seems 
that one could continue the process without limit. Going 
one or two more steps further, which we have not yet 
done, might be helpful for identifying a closed form ex- 
pression that would be positive to all orders. One might, 
after considerable effort, in this way find many expres- 
sions that had positive quasi-local energy for a region of 



any size. 

While we have found an infinite number of new quasi- 
local holonomic pscudotensor energy-momentum expres- 
sions with a certain desirable property, yet, as far as 
we have investigated them, none of these holonomic ex- 
pressions appears to be natural or especially appealing. 
We can offer an alternative: the teleparallel gauge cur- 
rent (see, e.g., 3) is naturally proportional to the Bel- 
Robinson tensor in the small vacuum region limit [32] . 

It should be remarked that our analysis also applies to 
many quasi-local gravitational energy expressions which, 
while not formulated within the pseudotensor framework, 
nevertheless take the same form in these limits. We have 
shown that, while the small region quasi-local energy pos- 
itivity Bel-Robinson requirement is a strong condition (in 
particular, it excludes all the classical pseudotensors), 
nevertheless it still allows for a lot of freedom in the 
choice of quasi-local energy-momentum expression. 
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